GW+?7 real-space Green's function calculations of x-ray spectra 

Towfiq Ahmed, J. J. Kas, and J. J. Rehr 

Dept. of Physics, Univ. of Washington Seattle, WA 98195 

(Dated: July 12, 2011) 

Abstract 

^^ . The Hubbard model is implemented in real-space Green's function calculations of x-ray spectra 

using an effective self-energy adapted from the LSDA+C/ method of Anisimov et al. This self- 
energy consists of an energy-dependent many-pole approximation to the GW self-energy with an 
additive correction due to on-site Coulomb repulsion among the partially filled localized-electron 
V , states. This leads to a GW+U approach which provides an efficient procedure to account for local 

H ' correlation effects on x-ray spectra. Results are presented for the spin and angular momentum 

projected density of states of MnO, NiO, and La2-xSrx.Cu04 (LSCO), for the K-edge x-ray spectra 
^ . of O atoms in MnO and NiO, and the unoccupied electronic states and O K-edge spectra of undoped 

'^ I LSCO. The method is found to yield reasonable agreement with experiment. 
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I. INTRODUCTION 

Density functional theory (DFT) together with quasi-particle corrections has been re- 
markably successful in describing the electronic structure and band-gaps of weakly interact- 
ing systems. For such systems, these quasi-particle corrections are often well described in 
terms of Hedin's GW self-energyii^ (where G refers to the one-particle Green's function and 
W the screened Coulomb interaction). On the other hand, the GW approach is generally 
inadequate to describe the band gap and other electronic properties in materials with well 
localized "id or 4/ electrons.- On the other hand, the strong local Coulomb interactions 
in these systems can often be approximated using a Hubbard-model,- in which the on-site 
electron-electron repulsion is represented by spin- and orbital-occupancy dependent "Hub- 
bard parameters" U and J . Combining the local spin density approximation (LSDA) of DFT 
with the Hubbard model leads to the LSDA+f/ method. In practice, the Hubbard correction 
is added to the original Kohn-Sham LSDA Hamiltonian while an approximate mean-field 
term is subtracted to avoid double-counting.- In order to calculate the excited state prop- 
erties and x-ray spectra of correlated systems, it is desirable to go beyond LSDA+f/ and 
incorporate energy dependent self-energy effects in terms of Hedin GW self-energy. This is 
the approach adopted here, which we refer to as the GW+f/ method. A related approach 
has been proposed by Jiang et al, where an accurate GW self-energy is calculated starting 
with LSDA+f/. Besides providing an improved screening model, their approach correctly 
predicts the band-gap in several d and / electron systems.-*^ In another prescription, Bansil 
et al developed a self-consistent GW+f/ scheme based on the tight-binding approximation 
and a single-band Hubbard model.-i^ Their method is found to qualitatively explain several 
pre-edge spectral features in high Tc cuprates.-i^^ Here we add Hubbard correction terms 
to an approximate many-pole GW self-energyii in a single-step approach, although for- 
mally such corrections could be incorporated within a self-consistent GW framework.- The 
implementation of our approach within the real-space Green's function (RSGF) formalism 
simplifies the calculation compared to conventional LSDA/GW+f/ methods, and is one of 
the primary goals of this paper. With the aid of this extension, we investigate the effects of 
correlated rf-electron states on the angular momentum projected density of states (/DOS), 
and the excited state spectra including x-ray absorption spectra (XAS) and x-ray emission 
spectra (XES) of a number of materials. As in other codes, e.g., WIEN2K— and SPRKKR,— 



the Hubbard parameters are here treated as fitting parameters. 

Our RSGF/GW+f/ method is tested on several (i-electron systems including MnO, NiO 
and the undoped high T^ cuprate La2_a;Sra;Cu04 (LSCO). In these materials, the electronic 
structure and band gaps are strongly influenced by U, the charge-transfer energy A, and 
the one-electron band-width W. Related calculations for MnO and NiO have also been 
carried out using a combined GW/LDA+f/ approach by Jiang et al^ Treatments of Ti oxide 
compounds using LDA+[/ within the multiple scattering formalism have also been carried 
out by Kriiger.— We flnd that our approach yields reasonable agreement with bulk-sensitive 
probes such as XES and XAS which are used to measure band gaps between occupied and 
unoccupied states.— 

II. THEORY AND METHODOLOGY 

In this Section we describe our implementation of the GW+f/ method as an extension of 
the real-space Green's function (RSGF) multiple-scattering formalism.— li^ Our implemen- 
tation follows the strategy used in the FEFF9 code and thus permits calculations of both 
electronic structure and x-ray spectra that can account for local atomic-correlation effects. 
Hartree atomic units {e = % = m = 1) are implicit unless otherwise specifled. 

A. RSGF Method 

We begin with a brief outline of the RSGF formalism used in this work. In this approach 
physical quantities of interest are expressed in terms of the local quasi-particle Green's 
function G{t,v' ,E). For example, the physical quantity measured in XAS for photons of 
polarization e and energy u is the x-ray absorption coefficient /i(ci;), 
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li{uj) (X Im (0c| e- rG(r, r',c<; + i?c)e ■ r'|0c) , (1) 
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where Ec is the core electron energy and |</)c) is the core state wave function. The FEFF9 code 
also calculates closely related quantities such as the spin and angular momentum projected 
density of states (/DOS) p|" (-E) at site n, 

pt\E) = --Im J2 / GZdr, r, E) r' dr, (2) 

i-n " 



where Rn is the Norman radius around the n*'' atom,— which is analogous to the Wigner- 
Seitz radius of neutral spheres, and the factor 2 accounts for spin degeneracy. The coefficients 
G^'*^, characterize the expansion of the Green's function G(r, r', E) in spherical harmonics, 

G(r,r',i?) = Y^ Y,{r)GZ,>{ry,E)Yl,{r'), (3) 

L,L',a 

where L = (/, m) denotes both orbital and azimuthal quantum numbers. In these formulae, 
the quasi-particle Green's function for an excited electron at energy E is given formally 
(matrix-indices suppressed) by 

g{e) = [e-h-j:{e)]-\ (4) 

where H is the independent-particle Hamiltonian 

H = ^ + V, (5) 

and V is the Hartree-potential. For convenience in our calculations, the Hamiltonian is re- 
expressed in terms of a Kohn-Sham Hamiltonian H^^ = H + Vxc where V^c is a ground state 
exchange-correlation^ functional, and the self-energy is replaced by a modified self-energy 
S(£') — Vxc which is set to zero at the Fermi-energy E = Ep. In this work we use the von 
Barth-Hedin LSDA functional V^.c[n(r),m(r)],— where n(r) = n^ + U]^ is the total electron 
density and ?7i(r) = n^^ — n^ is the spin polarization density. In practice, it is useful to 
decompose the total Green's function G{E) as 

G{E) = G\E) + G'''{E), (6) 

where G'^{E) is the contribution from the central (absorbing) atom and G'^'^{E) is the 
scattering part. Full multiple scattering (FMS) calculations can be carried out by ma- 
trix inversion, i.e., with G = [1 — G'°r]~^G°, where G^ is the bare propagator and T is 
the scattering T-matrix, which are represented in an angular-momentum and site basis: 

G^ = <^nL,n'L'(-^)[l " ^n,n'] and T = t'^LSi^i'6m,m'Sn,n'Sa,a' ■ Finally, f^^ IS the siugle site 
scattering matrix, which is related to the single site phase shifts, i.e., 

t;;i = exp(^5:jsin(5:;j. 

Within the spherical muffin-tin approximation, G^{E) can be expanded in terms of the 
regular Ri{y,E) and irregular Hi{r,E) solutions of the single site Schrodinger equation.— 



In the FEFF code a typical calculation of the electronic structure (ground or excited state) 
starts with a self-consistent calculation of the electron density and Kohn-Sham potentials.— 
Once the self-consistent potential is obtained, the Green's function is constructed and used 
to calculate XAS and other quantities of interest. Of particular interest in this paper is the 
local spin-dependent density matrix for the n-th site 

n""' ,,= --[ ^ dE [ IniG""' , ,(r,r, E)d\ (7) 

nlm,nlm -rr / / nlm,nlm \ ' ' / ' \ / 

^ J J cell 

where the n denotes the cell defined by the Norman sphere centered about the n*^ atom, r, r' 
are relative to the center of the cell -R„, and a is the spin-index, and we explicitly designate 
the azimuthal quantum numbers m and m' . For a more detailed description of the multiple 



scattering RSGF method see Refs. [20|J21|. 



B. GW-^C/ Self-energy 

Quasi-particle effects are key to an accurate treatment of excited state spectra,— and 
hence we need a good approximation for the electron self-energy. Current approximations for 
the self-energy typically begin with Hedin's GW approximation (GWA),- which is formally 
given by 

S = iGW (8) 

where G is the one electron Green's function, W = e~^v, is the screened-Coulomb interac- 
tion, and V the bare-Coulomb interaction. The FEFF9 code uses several approximations 
for the self-energy with the aim of providing efficient calculations of the energy dependent 
shift and broadening of spectral features. The default, which is appropriate at high energies, 
is the Hedin-Lundqvist plasmon-pole model,-i2^ based on the electron gas and a single-pole 
approximation to the dielectric function. An extension which improves the self energy at low 
energies is a many-pole model, where the the dielectric function is represented as a weighted 
sum of poles matched to calculations of the loss function in the long wavelength limit.— Al- 
though these models significantly improve quasi-particle calculations of unoccupied states, 
they do not necessarily obtain accurate band-gap corrections. In our implementation of 
GW+f/, an energy, spin and orbital dependent total potential is constructed that incorpo- 
rates the GW plasmon-pole or many-pole self energy T,{E) and the Hubbard correction V^J^, 
with parameters chosen to obtain the correct gap. Although such a construction can be 



done using self-consistent methods,— here we use only a single-step calculation. Thus we 
define our total potential as 

y(r, E) = V^'^'^^ir) + S«^(E) + V^^- (9) 

The orbital and spin-dependent Hubbard contribution to the potential is calculated as de- 
scribed in the next section. 

C. LSDA-I-C/ formalism 

Our construction of Vl^^{E) is adapted from the LSDA+[/ approach of Anisimov et 
al— In their approach one starts with the total energy functional of the system and adds 
a Hubbard correction to account for the Coulomb interaction between localized, strongly 
correlated electrons. It is generally assumed^^ that a similar mean- field term should exist 
in LSDA or other DFT approaches which must be subtracted from the energy functional to 
avoid double counting, 

E^K(f),n"] = E^^^^K(f)] (10) 

where n'^{r) is the charge density, n'^ the density matrix, E^ the Hubbard interaction, and 
Edc the double counting term. The Hubbard term depends on the density matrix n'^i'^^i^,, 
and on-site Coulomb interactions between the localized electrons. 

As discussed by Albers et al.,— an ab initio determination of the Hubbard parameters 
is not straightforward, since they are sensitive to screening of the Coulomb interaction. 
Hence, the Hubbard terms are often regarded as fitting parameters while the density matrix 
is calculated from first-principles to construct the Hubbard potential. This is the approach 
adopted here. While it may be possible to go beyond this parametrization and calculate the 
Hubbard terms using approaches such as constrained-LDA— i^^, or constrained-RPA,— 12^, 
such estimates are beyond the scope of this paper. For systems where the localized electrons 
are atomic-like, the density matrix can be approximated^^ as 

This spherical approximation is reasonable for many systems including TMOs, and good 
agreement for the band gap is found when the non-sphericity of d-d interaction as well 
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as the off-diagonal terms of rimm' are ignored.- With these approximations, the number of 
parameters is reduced to only two, namely U and J representing the screened direct and 
exchange intra-atomic Coulomb interactions, respectively. 
The total energy functional can then be written as 



E = E^^^^ + i Y. U{nl,-n'^){n~j-n°) 

m,m' ,a 

+ \ H {U-J){nl,-n°){nl,,-n"). (12) 

Here the double counting term E^c is represented by n° where n° = rid/lO and n^ = J2ma 'Kn- 
Using V{r) = 5E/5n(j{f), a simplified expression for the total LSDA+f/ potential is finally 
obtained,— i.e., 

where 

VL = U $^(n J, - n,) + {U-J)J2 K™' - ^o). (14) 

Within the spherical approximation, we need only consider the diagonal elements nf^ 
of the density matrix defined in Eq. (jTj). In a single-step spin-dependent calculation using 
the von Barth-Hedin LSDA functional, we first obtain nf^. In this prescription, a prior 

'f I 

knowledge of spin polarization of i-th atom rrii = nl — ni is required. For Mn, Ni, and Cu 
we used ?Ti = 5, 2, and 1 correspondingly using Hund's multiplicity rule^"^ for free atoms 
which is often treated as good approximations for such systems. 

The occupancy of the spin-up and -down states within the d-orbitals are thus determined 
in this single-step LSDA approach. Our calculations of spin-orbital occupancies of Mn and 
Ni (i-states using this scheme are listed in Tables I and II. Thus we essentially start with a 
spin dependent ground state calculation and introduce spin and orbital dependence using 
Anisimov's prescription of Hubbard model. This LSDA+f/ prescription is found to provide 
good agreement between the theory and experiment for the XAS of the TM compounds 
investigated here, although the self-consistent LSDA-|-f/ treatment may be more desirable 
in other cases. 

Values for the U and J parameters are taken either from previous work-i22, or chosen to fit 
the experimental band gap. For MnO and NiO, we used U ~ 4.5 and 7.5 eV and J = 0.9 eV, 
which are reasonably close to those calculated or discussed by other authors.-i^^ Using Eqs. 



(9) and (14) we then correct our self-consistent quasi-particle (QP) potential and obtain a 
new potential V'^^~^'^{r, E) given by 

yow+^ir, E) = V;^^^^(r) + S^^(i?) + V.^E). (15) 

Then using the GW+f/ Haniiltonian above, the wave functions RL{r,E) and Hl{t,E) 
are recalculated as solutions of the Schrodinger equation inside the muffin-tin spheres with 
our Hubbard modified potential. The orbital dependent phase shifts 5i^{E) are obtained 
by matching to the free solutions (spherical Bessel functions) at the muffin-tin, and the 
scattering t-matrices are found. 

C = exp(z(jrjsin(<5rj. (16) 

Finally the multiple-scattering equations are resolved with these t-matrices yielding the the 
total Green's function G = G"^ + G'^'^, which now includes the Hubbard-[/ correction. 
With the addition of the state dependent Hubbard correction, the potential of Eq. (15) can 
correctly account for the well known discontinuity^'^^ in exact DFT exchange-correlation 
potentials. However, such a term is absent from the conventional LDA and GGA approaches, 
rendering them incapable of including such band-gap corrections. 

III. RESULTS AND DISCUSSION 

A. Transition Metal Oxides 

Transition metal oxides (TMOs) such as MnO and NiO are considered to be prototypes 
of strongly correlated Mott type insulators with localized and partially filled d-electrons at 
the metal sites. These TMOs have NaCl like crystal structures, (Cubic 0\ symmetry, and 
fmSm space group). Below their respective Neel temperatures, they all exhibit a rhombohe- 
dral distortion due to anti-ferromagnetic (AF) ordering, which is also known as exchange 
anisotropy.— We examined the effects of such crystal distortions but they had negligible 
influence on the spectral features of interest here. In the following subsections we present 
results for the total and angular momentum projected DOS of MnO and NiO for a few 
values of U . The exchange parameter J is typically much smaller than U and variations 
were found to be small over the transition metals; thus we used J = 0.9 eV for all cases.- 
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For both compounds, the O K-edge XAS and XES are also calculated and compared with 
experimental results. 

1. MnO 

In order to compare with room temperature experiment,— we have taken an undistorted 
MnO crystal with a = h = c = 4.4316 A and a = (3 = 90.624°.— In this paper, we do 
not consider periodic magnetic effects; however, the single site moments are implicitly taken 
into account in our GW+f/ implementation. In our FMS calculations for MnO, we used a 
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FIG. 1: (Color online) U dependence on total DOS of MnO with spin up (solid red) and spin 
down (daslied blue) for different values of U: (a) GW {U=0), (b) U = 2.0 eV, (c) U = 4.0 eV and 
(d) U = 6.0 eV; Vertical dashed line is at the Fermi energy. 

cluster of 250 atoms, which was adequate to converge the spectrum, and a smaller cluster of 
60 atoms for the self-consistent muffin-tin potentials. For this system we calculated the O K 
edge XES and XAS and the spin and angular momentum projected DOS about the Mn and 
O sites with and without GW+f/ corrections. Fig. 1 shows a comparison of our calculated 
total ground state spin-resolved DOS of MnO to that calculated with different values of U. 
While a calculation with a GW self-energy underestimates an insulating gap (dashed blue 
line in Fig. 1), a gap close to that observed in experiment is obtained using the Hubbard 
correction U = 4.6 eV. When applied to Mn rf-states, an upper Hubbard band appears 
at about +1.5 eV, as seen in Fig. 2(a). The O p-states (Fig. 2(b)) near Ep are strongly 
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FIG. 2: (Color online) angular momentum projected Z-DOS for Mn and O in MnO with GW(L'^ = 
0.0) and GW+U {U = 4.6 eV); Spin up an down DOS are above and below the horizontal axis 
correspondingly: (a) Mn d-DOS GW(dashed) and GW+C/ (solid), (b) similar color and line style 
correspondence for O p-DOS; Vertical dashed line is at the Fermi energy. 

hybridized with Mn (i-states (Fig. 2(a)); thus a gap is also seen in the O p-DOS. However, 
the O j»-states around 6-8 eV only hybridize with Mn s-p-states (not shown) and are not 
affected by the Hubbard correction. In Table (I) we present the spin-orbital occupancies 
of the localized Mn d-states and the corresponding Hubbard correction for U = 4.6 and 
J = 0.9 eV. 

TABLE I: Mn d-state parameters {U = 4.6; J = 0.9 eV) 

1 m n,^ nL nL C(eV) ^i(eV) 

2 1.02 0.90 0.12 -1.82 1.61 
2 ±1 0.97 0.84 0.13 -1.56 1.55 
2 ±2 0.95 0.83 0.12 -1.52 1.63 



Bulk sensitive XES and XAS for TM oxides often provide a good assessment of the band 
gap in insulators.— In Fig. 3 we compare our GW+U calculation of the O K-edge XAS 
and XES with experiment.— Fig. 3 shows the result of our spin resolved FMS calculation 
obtained with both Hubbard and GW self-energy corrections (b) compared to results with no 
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FIG. 3: (Color online) O K-edge XAS (black) and XES (red) in MnO. (a) QP FEFF calculation 
using GW plasmon-pole self-energy, (b) GW+U {U = 4.6), and (c) experiment.— The vertical 
dashed lines are a guide to the eye. 

Hubbard correction (a) and experiment (c). The XAS calculation was done in the presence 
of a screened core-hole at the absorbing O atom while for XES no core-hole was included; 
these approximations are consistent with the final-state- and initial-state rules for XAS 
and XES respectively. Our Hubbard corrected self-energy blue shifts the first excitation at 
around 534 eV, while the rest of the unoccupied states, including the main peak at 540 eV, 
are unchanged. In XES, the highest occupied state moves down by 3 eV which is now on 
the other side of the second vertical dashed line in Fig. 3. These distinct, opposite shifts of 
the highest occupied and first unoccupied states are due to the strong hybridization of O 
p-states with the localized Mn d states. This can also be identified in Fig. 2(b) as the lower 
and upper Hubbard bands (UHB) at around -2 and 2 eV. 



2. NiO 



In order to compare with room-temperature experiments^ we have accounted for the 
rhombohedral distortion along the [111] direction.— i^ Our methods for calculating electronic 
structures of NiO are similar to those for MnO, except for the input NiO crystal structure, 
where we have used a slightly distorted crystal with a = b = 4.168 A, c = 4.166 A, and 
a = /3 = 90.055°, 7 = 90.082°. With the Hubbard correction, the best agreement with the 
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FIG. 4: (Color online) U dependence on total DOS of NiO with spin up (solid red) and spin down 
(dashed blue) for different values of U: (a) GW (C/=0), (b) U = 3.5 eV, (c) U = 7.5 eV and (d) 
U = 12.5 eV; Vertical dashed line is at the Fermi energy. 

experimental XAS was again obtained with U = 7.5 eV. Fig. 4 and Fig. 5 show the gap 
opening in the spin projected total DOS of NiO for higher values of U. The O p-states 
in NiO are also strongly hybridized with localized Ni rf-states as in MnO. The spin-orbital 
occupancies and corresponding Hubbard potential for the Ni d-states are listed in Table II. 

TABLE II: Ni d-state parameters {U = 7.5 eV; J = 0.9) 
^ ^ r^im nL nl V^|^(eV) l^t.(eV) 

2 1.22 0.97 0.25 -2.75 3.74 
2 ±1 1.64 0.96 0.68 -2.53 1.88 
2 ±2 1.44 0.97 0.47 -2.63 2.48 

Our GW plasmon-pole calculation in Fig. 6(a) exhibits considerable overlap between 
the O K-edge XAS and XES spectra, due to the underestimated insulating gap. However, 
the introduction of the Hubbard interaction {U = 7.5 eV) increases the gap, causing the 
pre-peaks of both the XAS and XES to split further apart, as shown in Fig. 6(b). For 
comparison, we also show a WIEN2K LDA+U calculation in Fig. 6(c) for the O K-edge 
EELS in NiO.^ 
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FIG. 6: (Color online) NiO O K-edge XAS (black) and XES (red) experiment vs theory: (a) FEFF 
GW+U plasmon-pole (PP) self-energy; (b) FEFF GW+U with Hubbard correction U = 7.5 eV; (c) 
GW+U calculatioi^ for O Is EELS using WIEN2K; and (d) Experiment .-IS The vertical dashed 
lines are a guide to the eye. 




E-Ep (eV) 

FIG. 5: (Color online) angular momentum projected /-DOS for Ni and O in NiO with GW{U = 
0.0) and GW+U (U = 7.5 eV); Spin up an down DOS are above and below the horizontal axis 
correspondingly: (a) Ni d-DOS GW(dashed) and GW+U (solid), (b) similar color and line style 
correspondence for O p-DOS; Vertical dashed line is at the Fermi energy. 

Aligning the first peak of this calculation with experiment [Fig. 6(d)], we observe an un- 
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derestimation of the high energy peaks at around 544 eV. These peaks can be attributed to 
O p-states which are strongly hybridized with Ni s— and p-states. Similar behavior has been 
found in NiO,-i^ and other TM compounds.— i^ We attempted to improve these results by 
using a more accurate GW many-pole self-energyii for NiO, while applying the Hubbard 
correction to the Ni (i-states. This many-pole self-energy includes a more realistic treatment 
of inelastic losses than the plasmon pole model, and yields improved agreement with exper- 
iment, as seen from Fig. 6(b). These results demonstrate that an accurate treatment of the 
delocalized s-p-states can also be important in such systems. Thus in order to achieve good 
agreement between theoretical and experimental spectral features, a systematic considera- 
tion of excited state properties including both locahzed- and delocalized states seems to be 
important. 

B. LSCO 

In recent years, understanding the doping dependence of high T^ cuprates has become an 
interesting challenge. LSCO (La2-xSr^.Cu04), which is a prototype of hole- doped cuprates, 
exhibits metallic and paramagnetic behavior at high doping,— and becomes an AF insulator 
when undoped. Between these limits, the system goes through a superconducting phase 
at about x = 0.15. A good description of the electronic structure in its insulating phase is 
important to understand the exotic doping dependent phase transformations in such systems. 

In the over-doped region with doping concentrations x > 0.2, LSCO becomes param- 
agnetic, and is well described by a self-energy approximation constructed from a single 
band Hubbard model.— A Fermi-liquid description thus becomes more appropriate for such 
systems. As doping is reduced, correlation effects due to localized states become more 
important, and the implementation of Hubbard U to the d electrons on the Cu sites is seen 
to open a gap. A gap correction using GW+f/ on the partial d-DOS of Cu and p-DOS 
of O is shown in Fig. 7. Our O K-edge XAS for GW and GW+f/ with U = 7.0 eV are 
compared with experimental results in Fig. 9. Our result agrees qualitatively with the 
undoped LSCO experiment, while the over-doped LSCO system is adequately reproduced 
by a GW calculation alone (f/ = 0). This result is not surprising, since in the absence of 
the Hubbard term, the LDA does not predict a correlation gap. As a result the system 
is predicted to be metallic, mimicking the over-doped (x ~ 0.3) paramagnetic phase of 
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FIG. 7: (Color online) (a) GW+U {U=7.5 eV) Cu spin up (solid red) and down (dashed blue) 
d-DOS; (b) GW (f/=0.0) Cu spin up and down d-DOS. The vertical dashed line is at the Fermi 
energy. 
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FIG. 8: (Color online) (a) GW+U (C/=7.0 eV) O spin up (solid red) and down (dashed blue) 
p-DOS; (b) GW (C/=0.0) O spin up and down |)-DOS. The vertical dashed line is at the Fermi 
energy. 

La1_a.Sra.CuO4. A complete description of the doping dependence of spectral features from 
over-doped {x = 0.3) to undoped (x = 0.0), requires a dynamical self-energy correction that 
incorporates pseudo-gap, superconducting, and Fermi-liquid physics.— 
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FIG. 9: (Color online) O K-edge XAS for LSCO: (a) our GW+U calculation with for U = 7.0 (red) 
and GW only (black); (b) experimental K-edge XAS for undoped (x = 0.0, red) and over-doped 
(x = 0.3, black) LSCO, and the vertical dashed lines are a guide to the eye. 

IV. SUMMARY AND CONCLUSIONS 

We have implemented a Hubbard model adapted from the LSDA+f/ method of Anisimov 
et al. as an extension of the real-space Green's function approach for calculations of x-ray 
spectra of correlated materials. In our construction two parameters U and J are chosen 
to match the experimental gap. These Hubbard parameters are introduced in terms of an 
effective self-energy correction leading to a GW+f/ approach which provides an efficient way 
to account for local correlation effects on x-ray spectra. Such a theoretical understanding 
of O K-edge XAS and XES is useful to explain key electronic features of strongly correlated 
systems For example, in the AF insulating phases of transition metal oxides, several impor- 
tant features of the experimental profile of the O K-edge XAS and XES, and in particular 
the correlation-gap, require theoretical treatments beyond the quasi-particle approximation. 

Our GW+f/ approach yields results which are in good agreement with experiment for 
the O K-edge spectrum of MnO and NiO. However, the agreement is only qualitative for 
more complex systems such as LSCO. This suggests the need for including a more com- 
prehensive treatment of superconducting and pseudo-gap physics capable of incorporating 
doping dependence in the under-doped regime of such systems.—"^ Finally we note that 
our current approach is limited to the quasi-particle approximation together with Hubbard 
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model corrections, while many-body effects such as satellites are neglected. However, some 
of these many-body aspects can be obtained by incorporating additional charge transfer 
contributions in the Hamiltonian.— i^ 

V. ACKNOWLEDGMENT 

We thank A. Bansil and R. Markiewicz and especially P. Rinke for stimulating sugges- 
tions. This work is supported by the Division of Materials Science & Engineering, Basic 
Energy Sciences, US Department of Energy Grants DE-FG03-97ER45623 and DE-FG02- 
07ER46352. This research also benefited from the collaboration supported by the Compu- 
tational Materials Science Network (CMSN) program of US DOE under grant DE-FG02- 
08ER46540. 



1 G. Onida, L. Reining, and A. Rubio, Rev. Mod. Phys. 74, 601 (2002). 

2 L. Hedin and S. Lundqvist, Solid State Phys. 23, 1 (1969). 

^ V. I. Anisimov, J. Zaanen, and O. K. Andersen, Phys. Rev. B 44, 943 (1991). 
^ V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, J. Phys.: Condens. Matter 9, 767 
(1997). 

5 H. Jiang, R. I. Gomez- Abal, P. Rinke, and M. Scheffler, Phys. Rev. B 82, 045108 (2010). 

6 H. Jiang, R. I. Gomez- Abal, P. Rinke, and M. Scheffler, Phys. Rev. Lett. 102, 126403 (2009). 
■^ T. Das, R. S. Markiewicz, and A.Bansil, Phys. Rev. B 81, 174504 (2010). 

8 R. S. Markiewicz, S. Sahrakorpi, and A. Bansil, Phys. Rev. B 76, 174514 (2007). 

^ S. Basak, T. Das, H. Lin, J. Nieminen, M. Lindroos, R. S. Markiewicz, and A. Bansil, Phys. 

Rev. B 80, 214520 (2009). 
^^ T. Ahmed, T. Das, J. J. Kas, H. Lin, R. S. Markiewicz, A. Bansil, J. J. Rehr, B. Barbiellini, 

and F. D. Vila (2011), in Press. 
^^ J. J. Kas, A. P. Sorini, M. P. Prange, L. W. Cambell, J. A. Soininen, and J. J. Rehr, Phys. 

Rev. B 76, 195116 (pages 10) (2007). 

12 O. Bengone, M. Alouani, P. Blochl, and J. Hugel, Phys. Rev. B 62, 16392 (2000). 

13 H. Ebert, A. Perlov, and S. Mankovsky Solid State Communications 127, 443 (2003). 



17 



^^ p. Kriiger, Journal of Physics: Conference Series 190, 012006 (2009). 

^^ E. Z. Kurmaev, R. G. Wilks, A. Moewes, L. D. Finkelstein, S. N. Shamin, and J. Kunes, Phys. 

Rev. B 77, 165127 (2008). 
^^ J. J. Rehr, J. J. Kas, F. D. Vila, M. P. Prange, and K. Jorissen, Phys. Chem. Chem. Phys. 12, 

5503 (2010). 
^'^ J. J. Rehr, J. J. Kas, M. P. Prange, A. P. Sorini, Y. Takimoto, and F. Vila, Comptes Rendus 

Physique 10, 548 (2009). 
^^ A. Ankudinov and J. Rehr, Phys. Rev. B 62, 2437 (2000). 

19 U. von Barth and L. Hedin, J. Phys. C: Solid State Phys. 5, 1629 (1972). 

20 A. Ankudinov, B. Ravel, J. Rehr, , and S. Conradson, Phys. Rev. B 58, 7565 (1998). 

21 J. J. Rehr and R. C. Albers, Rev. Mod. Phys. 72, 621 (2000). 

22 W. von der Linden and P. Horsch, Phys. Rev. B 37, 8351 (1988). 
2^ A. Georges, AIP Conference Proceedings 715, 3 (2004). 

24 R. C. Albers, N. E. Christensen, and A. Svane, J. Phys.: Condens. Matter 21, 343201 (2009). 
2^ O. Gunnarsson, O. K. Andersen, O. Jepsen, and J. Zaanen, Phys. Rev. B 39, 1708 (1989). 
26 I. V. Solovyev and M. Imada, Phys. Rev. B 71, 045103 (2005). 
2^ F. Aryasetiawan, K. Karlsson, O. Jepsen, and U. Schonberger, Phys. Rev. B 74, 125106 (2006). 

28 M. Cococcioni and S. de Gironcoh, Phys. Rev. B 71, 035105 (2005). 

29 V. I. Anisimov, I. V. Solovyev, M. A. Korotin, M. T. Czyzyk, and G. A. Sawatzky, Phys. Rev. 
B 48, 16929 (1993). 

^^ B. Fromme, Electronic Structure of MnO, CoO, and NiO, vol. 170 of Springer Tracts in Modern 

Physics (Springer Berlin Heidelberg, 2001). 
■^1 S. Liu and W. Langenaeker, Theoretical Chemistry Accounts: Theory, Computation, and Mod- 

ehng (Theoretica Chimica Acta) 110, 338 (2003). 
^2 F. Aryasetiawan, K. Karlsson, O. Jepsen, and U. Schonberger, Phys. Rev. B 74, 125106 (2006). 

33 J. P. Perdew, R. G. Parr, M. Levy, and J. L. Balduz, Phys. Rev. Lett. 49, 1691 (1982). 

34 A. E. Berkowitz and K. Takano, Journal of Magnetism and Magnetic Materials 200, 552 (1999). 

35 B. Morosin, Phys. Rev. B 1, 236 (1970). 

36 K. Nakahigashi, N. Fukuoka, and Y. Shimomura, J. Phys. Soc. Jpn. 38, 1634 (1975). 
3^^ H. Kondoh and T. Takeda, J. Phys. Soc. Jpn. 19, 2041 (1964). 

38 L. V. Dobysheva, P. L. Potapov, and D. Schryvers, Phys. Rev. B 69, 184404 (2004). 

18 



39 H. Kurata, E. Lefevre, C. Colliex, and R. Brydson, Phys. Rev. B 47, 13763 (1993). 

^° T. Kotani and M. van Schilfgaarde, J. Phys.: Condens. Matter 20, 295214 (2008). 

^1 R. S. Markiewicz, T. Das, and A. Bansil, Phys. Rev. B 82, 224501 (2010). 

^^2 T. Das, R. S. Markiewicz, and A. Bansil, Phys. Rev. B 81, 174504 (2010). 

^3 L. Hedin, J. Phys.: Condens. Matter 11, R489 (1999). 

''^ J. D. Lee, O. Gunnarsson, and L. Hedin, Phys. Rev. B 60, 8034 (1999). 



19 



